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ABSTRACT 



Context. Optical-infrared interferometry can provide direct geometrical measurements of the radii of Cepheids and/or reveal unknown 
binary companions of these stars. Such information is of great importance for a proper calibration of Period-Luminosity relations and 
for determining binary fraction among Cepheids. 

Aims. We observed the Cepheid X Sgr with VLTI/ AMBER in order to confirm or disprove the presence of the hypothesized binary 
companion and to directly measure the mean stellar radius, possibly detecting its variation along the pulsation cycle. 
Methods. From AMBER observations in MR mode we performed a binary model fitting on the closure phase and a limb-darkened 
model fitting on the visibility. 

Results. Our analysis indicates the presence of a point-like companion at a separation of 10.7 mas and 5.6 mag^ fainter than the 
primary, whose flux and position are sharply constrained by the data. The radius pulsation is not detected, whereas the average 
limb-darkened diameter results to be 1.48±0.08 mas, corresponding to 53±3 Rq. at a distance of 333.3 pc. 
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Classical Cepheids are fundamental astrophysical objects, be- 
ing the most popular primary distance indicators in the nearby 
Universe. These pulsating stars obey well defined Period- 
Luminosity (PL) relations, so that it is possible to derive their 
distance moduli by measuring periods and apparent magnitudes. 

Cepheids can be used to trace the evolutionary properties of 
Helium burning intermediate-mass stars, as they are crossing the 
instability strip along the so-called "blue loops" in the HR dia- 
gram. However, the comparison between pulsation and evolu- 
tionary masses disclosed the problem known as "Cepheid mass 
discrepancy" (ICoxlll980l) . due to the pulsation masses result- 
ing system a tically smaller than the evolu t ionary ones ( s ee e.g . 
Bono et all d200lh: iBeaulieu et al l (l200ll) : ICaputo et alJ <|2005); 
lKeller& Wood! (120061) 1. 

Recently two double eclipse binary Cep heids have been 
identif i ed in the Large Magel lanic Cloud dPietrzvnski et alJ 
d2011l) . iPietrzvriski et al.l d201 ll) ). providing the first measure- 
ments of the dynamical mass of classical Cepheids with 1-3% 
accuracy, whose values support the pulsation masses, or the 
Mass-Luminosity relations based on evolution ary models that 
account for mild convective core ov ershooting (ICassisi & Salarisl 
(1201 ll) : lPrada Moroni et al.l (1201 2h t. 

For non-eclipsing binaries, direct orbital determination of the 
companion by means of optical interferometry would provide 
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dynamical mass determinations able to disentangle the mass dis- 
crepancy problem. 

The binary fraction of the Cepheids has an impact on the 
Cepheid distance scale in the optical bands, since the compan- 
ions are typically main sequence (MS) stars. The binarity anal - 
ysis has been applied to a limited sample (Szabados 2003), 
but current estimates indicate that the bin ary fraction amon g 
MS B-type stars i s of the order of 60-70% dBrott et all (1201 ll) : 
IChini etaU (120121) 1. 

In order to improve the current observational scenario we 
have observed the Cepheid X Sgr with the VLTI. This partic- 
ular target was selected becau se: i) its trigonometri c parallax 
has been recently measured by Bened ict et al.l d2007l): i i ) accu - 
rate radial velocity curves are available (IMathias et alj |2 006): 
IStorm etaf] (1201 ll) ). and Hi) it was suggested bv ISzab ados 
(1990), on the basis of temporal changes in the y-velocity, that 
X Sgr might b e a binary, although t his finding has alway s been 
controversial dMathias et alJ d2006l) : lGroenewegenl (12008)). 



2. Observations and data reduction 

2.1. Data acquisition 

We observed X Sgr (HD 16 1592, R.A. 1 7:47:33.62, Dec. - 
27:49:50.83, J2000) using VLTI/ AMBER dPetrov et all 12001 
with the Auxiliary Telescopes. We used medium spectral resolu- 
tion in the jf-b and (MR, R~1500, 2.1pm) and the fringe tracker 
FINITO dGai etalj f2004) with a frame exposure time of 1 s. 
Observations were carried out with a maximum projected base- 
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Table 1. Log of the interferometric observations, with details on the SCIence-CALibrators sequence used and the optimal frame 
selection adopted for V 2 and (p. Magnitude and diameters of the two calibrators employed are shown. 



Date 


Target sequence 




Baselines 


% of frames for V 2 


% of frames for <f> 


2011-04-24 


CAL 2 -SCI-CAL 2 


-CALrSCI-CALi 


D0-I1-G1 


10% 


100% 


2011-05-03 


CAL 2 -SCI-CAL, 


-SCI-CAL r CAL 2 -SCI-CAL 2 


D0-A1-C1 


30% 


100% 


2011-05-07 


CALi-SCI-CAL] 


-SCI-CALi 


D0-A1-C1 


10% 


100% 


2011-06-11 


CALpSCI-CAL! 




K0-A1-I1 


50% 


100% 


2011-06-25 


CALrSCI-CAL, 


-SCI 


K0-A1-I1 


20% 


100% 


2012-03-31 


CAL 2 -SCI-CAL, 


-SCI-CALj 


K0-A1-G1 


20% 


100% 



ID 


Name 


mag£ 


9 a UD [mas] 


CALj 


HIP88839 


2.070 


1.881 ±0.133 


CAL 2 


HD157919 


3.203 


0.947 ± 0.066 



" Data from JMMC Stellar Diameters Catalogue - JSDC dLafrasse et al.ll2010h 



line length of 127.8 m, corresponding to 3.1 mas spatial resolu- 
tion, with the atmospheric seeing varied in the 0.7" -r 1.0" range. 

For the sake of clarity, we distinguish here a Run A (six 
nights from 201 1-04-24 to 201 1-06-25) and a Run B (2012-03- 
31 night only, see Table Q] and Fig. [TJ. For each run we used a 
CAL-SCI-CAL observation scheme (Table Q]) to better sample 
the variable instrument response during the night. 

2.2. Data reduction 

We u sed the amdlib v£] dTatulli et alj d2007l) . IChelli et alj 
(2009)) reduction package in conjunction with custom IDL 
scripts for optimizing the fringe selection, computing averages, 
and calibrating the visibilities. 

Each AMBER file contains many individual frames, which 
must be quality-selected to reject tracking losses in both FINITO 
and telescope guiding. Such events are associated to large val- 
ues in the optical path difference (OPD) and cause a poor sig- 
nal to noise ratio (SNR) of the fringes. We adopted the amdlib 
"PISTON+SNR" selection criteria, limiting the OPD variation 
below 8^m and keeping only the frames with the highest fringe 
SNR. The adopted percentage of frames, in the range 10 -r 50% 
for visibility V 2 and 100% for closure phase <f>, has been opti- 
mized for each night separately (see Table[TJ. 

Hence, we computed, for each file, the average values of V 2 
and </>, weighted by the error estimate provided by amdlib; the 
same was done for the calibrators. 

Finally, a night transfer function (TF) was computed using 
the calibrators listed in Table Q] 

3. Results and Discussion 

3. 1 . Visibilities and closure phases 

Fig. [1] shows the calibrated V 2 as a function of the spatial fre- 
quency B//1, while Fig.|2]displays the calibrated <p versus wave- 
length A. 

We see that the source is resolved, because V 2 is less than 
unity, and that significant deviations from a pure Uniform Disk 
(UD) are evident. This could indicate a possible detection of the 
Cepheid pulsation, or might suggest a more complex morphol- 
ogy, like that of a binary. In any case, in the <p plot we only 
see small deviations from zero, which is an indication that the 
companion, if present, should be faint and close to our detection 
limit. 

1 Available at http://www.jmmc.fr/amberdrs 



Both plots, in particular the V 2 , show that the error bars are 
clearly larger than the dispersion of the data points, indicating 
a correlated noise caused by the low-frequency variations (on 
time-scales of minutes) of the TF during the night. 

3.2. Model fitting 
3.2.1. Binary companion 

ISzabadosi (Il990h suggested, on the basis of changes in the y- 
velocity, that X Sgr might be a binary system with a period 
of ~507 days. The complex nature of X Sgr was further sup- 
ported by Sass elov & Lester! (Il990h who spectrally found that 
line asymmetries and line-doubling were well beyond the typi- 
cal Cepheids behaviour; in particular, the two absorption com- 
ponents showed different Doppler shifts, but they did not dis- 
play secular velocity changes compared with the y- velocity. This 
complex dynamical structure was also supp orted by other high - 
resolution (R~ 120,000) optical spectra bv lMathias et alj d2006). 
They hypothesized a binary, or even triple, nature of X Sgr. 
The d ata availa ble in the literature were analyzed bv lFeast et al.l 
(2008) and bv lGroenewegenl d2008) who also found evidence of 
binarity, these latter suggesting a slightly longer orbital period 
(-570 days). 

In order to check the companion hypothesis, we performed 
a binary model fit to the closure phase </>. We decided to use 
only this observable, for its low sensitivity on the unknown ra- 
dius variation, for which we adopted a crude UD estimate of 
fuD = 0.74 mas, provided by the V 2 values at the largest base- 
lines, namely the data from nights 201 1-04-24, 201 1-05-03 and 
2012-03-31. 

We started with Run A alone assuming no temporal variation 
in the binary configuration within the two months, which is a 
plausible assumption, considering the mentioned orbital period 
estimates. Therefore we fit the Run A data all together in order 
to increase accuracy. Run B was not included here, since it was 
acquired 10 months later. 

We adopted a grid search on a wide parameter space in 
Ax, Ay,f and r%, where Ax, Ay are the R.A. and Dec. offsets of 
the companion (these ar e searched within a radius of 100 mas 
around the primary, see lAbsilet al.1 (1201 hY ). / is the flux ratio 
of the companion to the primary, and is the companion UD 
radius. 

We obtain the best fit result for a flux ratio of / = (6± 1)- 10 3 
(corresponding to a 5.6 + 0.2 magnitude difference), while ri 
tends to zero, thus confirming the presence of a faint and un- 
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Fig. 1. Left: calibrated V 2 versus spatial frequency B/A for the entire set of observations, superimposed to the average UD model 
(continuous line). Right: spatial frequencies («,v-plane) coverage. Different colors and symbols refer to different observing nights 
(see legends). Note: the 201 1-06-25 night is not plotted in the left panel, due to a very unstable TF for V 2 . 
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Fig. 2. Calibrated </> versus A for the entire set of observations. Colors and symbols are the same as in Fig.Q] The yellow line shows 
the binary model described in Sec. l3.2| providing the best fit to the Run A observations. The dates reported above the panels refer to 
the observing nights. 



resolved companion, whose position is constrained by a single 
and well defined x 2 minimum at 7.15+0.04 S, 8.00+0.02 E mas 
from the primary (see Fig [3j. The parameters of the companion 
are summarized in Table [2] 

The fact that the best fit reduced x 2 is significantly smaller 
than one is due to the large error bars associated to the TF, al- 
ready described in Sec. 13.11 

To double check our finding, we built a series of synthetic 
binary models with different flux ratios and the same error bars 
of the real data, from which we derived a sensitivity limit of 
Slim — 2 • 10~ 3 , i.e. 3 times lower than the measured flux ratio. 

To refine this result we removed our assumption of fixed bi- 
nary configuration, in order to check if we could detect possible 
displacements among different nights. 

Therefore, we repeated the <p fitting separately for each night 
of Run A, to find out if the location of the y 1 minimum moves 
as a function of time (nights 201 1-05-03/201 1-05-07 and 201 1- 
06-1 1/201 1-06-25 were fitted as single points). 

The corresponding x 2 maps (Fig. [3] panel A) show several 
local minima within the 99% confidence interval, all of them 
having the same probability to be the true position of the com- 
panion. This shows us that a single baseline configuration is not 
sufficient to unambiguously locate such a faint object. 



However, each x 2 map shows a particular local minimum in 
the same position found using the entire data set, and even more 
importantly the little displacements of these minima are in the 
correct temporal sequence, as disclosed by Fig. [3] panel B. 

This empirical evidence further confirms our detection of the 
postulated companion of X Sgr, and although current data do not 
allow us to constrain the orbit, we have a hint on the direction of 
its motion. 

The binary fit for the single night of Run B also show the 
same kind of multiple minima x 2 map (not displayed), but with 
no minima close to the quoted position, as expected after a time 
interval of ten months. 



3.3. Angular diameter 

To determine the mean Cepheid diameter, and possibly its vari- 
ation, we fit the V 2 separately for each night by fixing the bi- 
nary parameters (Table |2]i and varying only the primary ra- 
dius r\. Following iKervella et alj d2004h. we ado pted the Limb 
Darkening (LD) intensity profile by Claret (2000). 

The radius r\ only depends on the mean value of each visi- 
bility spectrum (< V 2 (A) >), whose error bar we computed as the 
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Fig. 3. Portion of the reduced^ 2 map on the R.A. and Dec. offsets obtained by A) fitting the (p of the Run A data sets together with 
a point-like companion model at fixed flux ratio / = 6 • 10~ 3 , and B) for individual observing nights: red for 201 1-04-24, green for 
201 1-05-03/07, and blue for 201 1-06-1 1/25. The 1 4- 3cr confidence contours colours range from yellow to orange. 



Table 2. Parameters of X Sgr companion, and LD diameters. 
Random errors (first ones) and systematics (second ones) are 
distinguished. 



Ax [mas] 


Ay [mas] 


r 2 [mas] 


/ 


8.00+0.02 E 


7.15±0.04 S 


<0.5 


(6±1)-10- J 


Phase" 


Date 


9 LD [mas] 




0.83 


2011-05-03 


1.47+0.01+0.04 




0.54 


2011-04-24 


1.36±0.09±0.04 




0.35 


2012-03-31* 


1.59±0.08±0.04 





" Th e pulsation phas e was estimated using the formula provided 
by ISzabadoa d 19891) for constraining the pulsation period and the 
epoch of maximum of the different data sets. 

* For the 2012-03-31 night a pure LD model is adopted, given the 
undetermined binary parameters. 



quadratic sum of the weighted average error plus the correlated 
component due to the TF variatiorQ. 

For the diameter determination we only took into account 
the nights 201 1-04-24, 201 1-05-03 and 2012-03-31, which have 
sufficiently long baselines and good quality. The night of 201 1- 
05-03 plays a key role in this context, since it shows a pure ran- 
dom error with a negligible correlated component. 

Table [2] provides, in the bottom section, the best fit LD di- 
ameters a s a function of th e Cepheid pulsation phase, computed 
following lSzabadosi d 19891) . 

In the table we distinguish the random errors and the sys- 
tematics coming from the uncertainties in the calibrator diame- 
ters (TableQ}. Note that the latter are larger than the former only 
for our best data set, namely the night 201 1-05-03, showing the 
importance of using many different calibrators when a very ac- 
curate diameter measurement is required. 

The differences in the three diameter estimates, within ~ 1 <x 
from each other, do not allow us to positively detect the radius 
pulsation, and are probably due to the mentione d TF temporal 
variati on. The same limitation was also found by iKervella et alj 
(2004) using interferometric data collected with VINCI at VLTI. 

2 where the latter was estimated, for each visibility spectrum, as the 
squared difference between the r.m.s. of the error bars and the r.m.s. of 
the residuals with respect to a smooth polynomial fitting. 



Table 3. Recent X Sgr radius estimates. 



Article 


Radius [R Q ] 


This paper 


53 


± 3 


This paper 3 


55 


± 4 


Kervella et al. (,2004) 


53 


± 3 


Storm etal. (201 1? 


50 


± 3 


Groenewegen* c 


49 


± 3 



" Solution obtained by neglecting the presence of the companion. 

c IRSB method. 

c Private communication. 



Thus we only computed a mean physical radius of 53+3 R Q 
by adopti ng the trigonometri c parallax of 3.0+0.18 mas from 
the HST dBenedict et al.ll2007h . In Table [3] we co mpared it with 
recen t estimates from both optical interferometry (IKervella et alj 
2004) and Infra Red Surface Brightness (IRSB, Groenewegen 
(personal communication) and lStorm et al.l (1201 ll) ~). having used 
the same distance for all of these measurements. Finally, we also 
report the radius obtained by neglecting the binary companion, 
which has a ~3% impact on the mean radius. 



4. Conclusions 

We found evidence of the binary nature of the Cepheid X Sgr 
by analysing new interferometric data collected with AMBER at 
VLTI. The main conclusions of our work can be summarized as 
follows: 

i) Our analysis of the closure phase indicates that both the 
flux and the position of the companion can be constrained by 
our data sets. We find that the angular separation is 10.7 mas and 
that the companion is 5.6 mag fainter than the primary in the 
K-band. 

ii) The data allow us to directly measure an average Limb- 
Darkening diameter of 1.48+0.08 mas (i.e. a physical radius of 
53+3 R Q ) which agrees, at lcr level, with similar interferometric 
and IRSB mean radii for X Sgr available in the literature. 

Hi) We find that MR AMBER closure phases obtained from 
3 baseline configurations have the sensitivity needed to detect 
companions around nearby classical Cepheids. This evidence 
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becomes even more compelling if we take into account that cur- 
rent data were collected in medium quality seeing conditions. 

iv) The average visibilities of the 2001-05-03 night show 
that a non systematic error as small as 1 % can be obtained with 
AMBER in MR mode in stable TF conditions. This would be 
sufficient to clearly trace the differential diameter pulsation in 
future measurements, provided that the same calibrator star is 
used for all the observations. 

v) We find that, in the case of a binary Cepheid, the effect in- 
troduced by the companion on the primary radius measurement 
(~3%) is comparable with the expected radius variation (~5 -f- 
10%). 
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